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Abstract. We argue that the idea that the parton system created in relativistic 
heavy-ion collisions is formed in a state with transverse momenta close to 
thermodynamic equilibrium and its subsequent dynamics at early times is 
dominated by pure transverse hydrodynamics of the perfect fluid is compatible 
with the data collected at RHIC. This scenario of early parton dynamics may help 
to solve the problem of early equilibration. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

It is now commonly accepted that the evolution of matter created in heavy- 
ion collisions at RHIC energies is most successfully described by hydrodynamics of 
the perfect fluid [TJ [H [3J. In particular, the hydrodynamical approach reproduces 
reasonably well the particle transverse-momentum spectra and the elliptic flow 
coefficient V2- Very recently, it has been also shown that the hydrodynamic model 
can describe consistently the HBT radii [4|. 

On the other hand, the hydrodynamic picture is challenged by a serious problem. 
To reproduce correctly the data, the hydrodynamic evolution should start at the 
time well below 1 fm after the collision takes place. Such fast equilibration time 
is not easy to achieve with elastic perturbative cross-sections, hence the success of 
the hydrodynamic approach inevitably leads to the puzzle of early thermalization. 
This puzzle was widely discussed and several exotic mechanisms were proposed for its 
solution [51 [S] , however no one was yet accepted as fully satisfactory. 

In the present paper we analyze the possibility that, at the early stages of high- 
energy collisions, the hydrodynamic evolution applies only to transverse degrees of 
freedom of the partonic system. The longitudinal dynamics is essentially described by 
partonic free-streaming (we consider here the midrapidity region, y rs 0). 

The main reason for our study is the observation that the early equilibration, if 
at all possible, is particularly difficult to achieve in longitudinal direction. On the 
other hand, the equilibration of the partonic transverse-momentum spectrum may 
be obtained much easier. In fact, it is well known that the transverse momentum 
spectra observed in nucleon-nucleon collisions are well described by the Boltzmann 
distribution [8]. 
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Figure 1. Transverse-momentum spectra of 7r + measured by PHENIX in 
the centrality class 30-40% (solid line) |17| and the model spectra of gluons 
for various choices of Ti and Tf . The lowest dash dotted curve corresponds 
to Ti = 200 MeV and T f = 180 MeV. The highest dash dotted curve was 
obtained for T = 300 MeV and 2> = 180 MeV. The four almost parallel 
lines represent our results for T = 250 MeV and for four final values of 
the temperature: Tf = 200 (long dashed line), 180 (dash dotted line), 
160 (dashed line) and 140 (dotted line) MeV. Note that the original n + 
experimental spectra collected at -^/sjvjv = 200 GeV have been multiplied 
by a factor of 3 to account for the total hadron multiplicity. 



In our approach the partonic system may be treated as a superposition of 
independent transverse clusters. The clusters are formed by the particles having the 
same value of the rapidity. The space-time evolution of the system is described by the 
hydrodynamic equations introduced in Ref. [9]. They are based on the energy and 
momentum conservation laws 

d^ v = 0, (1) 

where the energy-momentum tensor is defined by the formula 

T"" = n °" gT (WIT - - V"V) . (2) 

27TT 

In Eq. ([2]) the parameter no is the density of clusters in rapidity space, v g is the 
degeneracy factor (we assume that our system is dominated by gluons, hence we use 
v g = 16), T is the temperature, and r = \/t 2 — z 2 . The four vectors and V M are 
defined by the equations: 

= (uo coshf?, u x , u y , uq sinh^), (3) 

V = (sinhr?, 0,0, cosh rj), (4) 

where r\ is the spacetime rapidity, and u°,u x , u y are the components of the four- velocity 
of the fluid in the rest-frame of a cluster, where y = = 0, 



u 



31,0) = (u°,u x ,uy,0). (5) 
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Figure 2. The elliptic flow coefficient V2 as the function of the transverse 
momentum. The PHENIX experimental results for pions and kaons in the 
centrality class 20-40% and for the collision energy ^/snn = 200 GeV [18] 
are compared to the model calculations with T; = 250 MeV and for four 
final values of the temperature T f = 200, 180, 160 and 140 MeV. The curves 
are denoted in the same way as in Fig. 1. The best agreement is obtained 
for Tf m 180 MeV. 



One of the interesting features of our approach defined by Eqs. |1| - is that it 
naturally leads to the entropy conservation, which takes the form 

= 0, (6) 

where 

S» = — — — — ^— — ■ . (7) 

2ttt v ' 

Eq. (|6|) follows directly from the projection of Eq. {T]) with the energy-momentum 
tensor ([2]) on the four- velocity U v . More information concerning the formal aspects of 
transverse hydrodynamics may be found in Refs. [TD|[TT]. Here we only note that our 
approach differs from the similar model by Heinz and Wong [12] . where the entropy 
is not conserved and the system is treated as dissipative. 

To solve numerically the hydrodynamic equations |T]) we have to specify the initial 
conditions, appropriate for the physical situation encountered in gold on gold collisions 
at RHIC energy. In our approach we assume that the transverse profile of the initial 
energy density is given by the density of wounded nucleons [T3] (see also [21 [T5]V 
This density is determined, for a given centrality, from the Glauber model. In the 
theoretical calculations we consider the centrality class 20 — 40%, corresponding to 
the impact parameter b ~ 7.9 fm. The initial temperature Ti at the origin, x±_ — 0, is 
taken as a free parameter. 

The parton spectra and the elliptic flow coefficient v-i are evaluated using the 
Cooper-Frye prescription [16] 

dN uq v g 



d 2 P± dy (2tt) 2 



dZ^xW F(x,p), (8) 
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where £ is the hypersurface where the purely transverse evolution comes to the end, 
and F(x,p) is the underlying phase-space distribution function, see Ref. [S . In 
our calculations £ is determined by the condition of constant temperature T — Tf. 
The details of our calculation of the parton transverse-momenta and the elliptic flow 
coefficient v-i are given in Ref. |10| . 

Our results obtained for several choices of the parameters T and Tf are shown in 
Figs. 1 and 2. We find that with T ~ 250 MeV the slope of the transverse momentum 
spectrum may be easily adjusted to the one measured for the pion spectra [PT] . On 
the other hand, the experimental values for V2 [H| select the final-temperature value 
Tf = 180 MeV. This value of the final temperature is reached by a system in a 
relatively short evolution time of about 4 fm. The short evolution time implies that 
the transverse size of the system is relatively small, smaller than that obtained from 
the HBT measurements, thus consistently leaving space for further, three-dimensional 
expansion of the system. 

We thus conclude that the concept of the initial transverse equilibration and 
hydrodynamic evolution of the partonic system is compatible with the data and the 
proposed model represents a possible solution the the problem of early equilibration. 

Very recently, the successful description of the RHIC data has been achieved 
in a model which assumes initial free-streaming in three space directions followed 
by a sudden transition to standard hydrodynamics [I]. Similarly to our arguments, 
that model indicates that the successful description of the data may be achieved 
in the approach without complete thermalization in the early stage. Thus, further 
investigations of early non-equilibrium dynamics (for example, see Ref. |19] ) are 
certainly necessary and interesting - they bring us closer to the solution of the early 
thermalization puzzle. 
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